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Abstract 



Muon pair production in the process e+e~ e+e~/i+/i~ is studied using the 
data taken at LEPl {^/s ^ rriz) with the DELPHI detector during the years 
1992-1995. The corresponding integrated luminosity is 138.5 pb~^. The QED 
predictions have been tested over the whole range accessible at LEPl (from 
several GeV^/c^ to several hundred GeV^/c^) by comparing experimental distri- 
butions with distributions resulting from Monte Carlo simulations using various 
generators. Selected events are used to extract the leptonic photon structure 
function F2 . Azimuthal correlations are used to obtain information on addi- 
tional structure functions, and Fg, which originate from interference terms 
of the scattering amplitudes. The measured ratios Fj/ F2 and F^/F2 are sig- 
nificantly different from zero and consistent with QED predictions. 



(Eur. Phys. J. C19(2001)15) 
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1 Introduction 

The study of the process e~^e~ e^e~ n~ provides a good way to test QED up to 
the fourth order of a. The photon structure can be studied by measuring photon structure 
functions which can be extracted in the so-called "single tagged" mode (Fig. 1), where 
one of the scattered electron^ is detected ("tagged") in an electromagnetic calorimeter 
while the other scattered electron goes undetected ("untagged"). This process can also 
be used as a reference one for studies of the hadronic structure function of the photon, 
providing a basis for a better understanding of the detector performance and for checking 
the analysis procedure. 

Previous measurements of muon pair production in both the single tagged mode and 
the double tagged mode (where the scattered electron and positron are both detected) 
have shown good agreement with QED predictions [|l],|,|,|J|J|J^,|,|,|10],|lT],||,|l|,||], with 



one exception where an excess of data events was observed in the double tag mode. 

This study, based on the data collected by the DELPHI collaboration at LEPl at 
centre-of-mass energies from 89.4 to 93 GeV, complements those results. It improves on 
previous DELPHI measurements of the leptonic photon structure function F2 by in- 
cluding all the LEPl statistics and increasing the coverage by an order of magnitude. 
This paper also presents results of studies of the azimuthal correlations, which are used 
to extract the ratios of the structure functions F]^/ F2 and F^/ F2. 



2 Event kinematics 

In the single tagged mode, where the tagged and untagged electrons are scattered with 
polar angles Otag and Ountag and energies Etag and Euntag respectively, and the probe and 
target photons have four- momenta q = {xtagEbeam,q) and p = {xuntagEbeam,p), the cross 
section of the reaction e~^e~ — > e~^e~X is given by fTT 



d^a _ dn{xuntag) d'^a{e'y eX) 



dxdQ dXiiniag dXiintag dxdC^ 

These two factors, the flux of target photons and the 67 cross section respectively, are 
given by: 

dn(^Xuntag) 01 I /i , /-i \2\ i ( '^EbeamiX ^untag) . ^untaq\ , \ 

- ' 1 + (1 - Xuntag) ) m Sm — " 1 + ^untag , 



^^untag ^^untag ( \ '^e^untag 

(2) 

dMn - eA-) _ 2.c^ { (1 + (1 _ ,)^) f p.) _ ,^f2(:r. Q^ P'')} . (3) 



dxdQ"^ xQ 

Here F2 and F^ are structure functions of the photon, a is the QED coupling constant, 
= —q^ ~ AEtagEbeamS^^'^i^^tag/'^) IS the squarcd 4-momentum transfer, = — is 
the virtuality of the target photon, and x and y are the Bjorken variables 

^ = TT— - = , , y = -T.- — = 1 COS — (4) 



2q-p W^^ + Q^ + P^' ^ p-k SXXuntag Ei 



beam 



where = (q' + pY is the invariant mass of the 77 (or yU'^/x ) system, k is the initial 
four- vector of the tagged electron, and s = AEl^^^. 

^Throughout the paper the term "electron" will be used for the tagged electron/positron 
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The structure function F2 can be extracted from the dependence of the cross section 
on X and . But Fl is small and is weighted by the small factor y^, making its direct 
measurement impractical. 

However, additional structure functions can be studied by looking at azimuthal cor- 
relations of the final state particles. The differential cross section of the process can be 
written as 

d^a{e^ e/i+/i") 2Txa^ 1 + (1 — yY 
dxdyd cos 6* dx/4:TT Q'^ xy 

|(2xFt + eiy)FL) - p{y)FACOsx+ \t{y)FB cos2x^ , (5) 

where x is the azimuthal angle, defined in the 77* centre-of-mass frame as the angle 
between the planes formed by the photon axis and the muon and the scattered electron 
respectively (Fig. 2), and 6* is the angle between the muon and the photon axis. The 
functions p{y) and e{y) are given by p{y) = (2 — y)\/l — y/{\ + (1 — yY) and e{y) = 
2(1 — ?/)/(! + (1 — yY) fl^ and can be taken equal_to 1 in the accessible kinematical 
region. The differential structure functions Ft, Fl, Fa, and Fb give the corresponding 
standard structure functions Ft, Fl, Fa, and Fb after integrating appropriately over 
cos^* (see section 7.2) taking into account that Fa is antisymmetric in cos^* [0. The 
cross section can then be written as 

d^a{e'y e/i+^u") ^ 27ra^ 1 + (1 - yf 
dxdydx/2n xy 

n{^- {Fim) cos X + ^ {nin) cos 2x) . (6) 

The structure functions F^ are combinations of transition amplitudes for the different 
helicity states of the photons. The structure function is related to the interference 
term between the two transverse helicity states of the photons. It is identical to F2, 
which is related to the longitudinal polarization of the virtual photon, in leading order 
and for massless muons. 



3 DELPHI detector 



The DELPHI detector has been described in detail elsewhere [p!9| , p0| . In this analysis, 
the scattered electron was tagged using 

• the Small Angle Tagger (SAT), the main luminosity monitor during 1991-93, covering 
polar angles from 2.5° to 8° (172° to 177.5°); it was made of alternating layers of 
lead sheets (0.9 mm thick) and plastic scintillator fibres (1 mm in diameter), aligned 
parallel to the beam; 

• the Small angle Tile Calorimeter (STIC), the main luminosity monitor since 1994, 
covering polar angles from 1.7° to 10.3° (169.7° to 178.3°); the STIC is a sampling 
calorimeter with 49 sandwiches of 3.4 mm steel-laminated lead plates and 3 mm 
thick scintillator tiles giving a total thickness of ~27 radiation lengths; 

• the Forward ElectroMagnetic Calorimeter (FEMC) covering from 10° to 36.5° 
(143.5° to 170°) in polar angle, consisting of two 5 m diameter disks containing 
a total of 9064 lead glass blocks. 



3 



The energy resolution of the tagging calorimeters was around 5% in SAT and FEMC and 
3% in STIC for an incident electron energy of 45 GeV. 

For muon identification, DELPHI contained barrel and forward muon detectors, each 
consisting of at least 4 layers of drift chambers. The muon chambers covered 78% of the 
solid angle. 

Combining the information from the tracking detectors, the relative momentum reso- 
lution Up/p varied from 0.001 xp to 0.01 xp {p in GeV/c), depending on the polar angle 
of the charged particle. 



4 Monte Carlo simulation 



Two event generators were used in order to simulate the signal process e+e^ 
e^e~ BDKRC [21| which includes only the multiperipheral diagram (Fig.l) to- 
gether with QED radiative corrections, and DIAG36 which lacks the QED radiative 
corrections but includes also the bremsstrahlung, annihilation and conversion diagrams. 
DIAG36 was used to check the role of these additional diagrams. 

Several generators were used to estimate the backgrounds to the process studied: 
BDKRC m was used to simulate e+e" e+e~r+r~, TWOGAM p3[ to simulate hadron 
production in two-photon collisions, DYMU3 |24] for the e^e~ 
KORALZ p5[ for e+e" ^ r+r-(7). 



yU+yU (7) process, and 



The generated events were passed through the full simulation of the DELPHI detector 
and reconstructed using the same program as for the data. 



5 Event selection and correction 

Events were selected as single tagged dimuon candidates if the following requirements 
were met. 

• There was a cluster in one of the electromagnetic calorimeters with an energy deposi- 
tion greater than 0.6 xi?;,eam (hereafter called the tagged electron). If the cluster lay 
within the polar angle range 20° - 160°, it was linked to a detected charged particle. 

• There were exactly two additional particles with opposite charges and polar angles 
between 20° and 160°. The relative errors on their momenta were less than 1. Their 
impact parameters with respect to the average interaction point were below 4 cm 
in the transverse plane and 10 cm along the beam. Their track lengths seen in the 
tracking detectors were at least 30 cm. Their momenta were above 0.5 GeV/c and 
2.5 GeV/c and the sum of their momenta was below 30 GeV/c. 

• At least one of the additional particles with a momentum greater than 2.5 GeV/c 
was identified as a muon by the DELPHI standard muon tagging algorithm [ pO[] . 

• The invariant mass of the two additional particles was above 1.7 GeV/c^. This 
requirement reduced the contribution from diagrams other than the multiperipheral 
one to below 0.25% for the low and 2% for the high sample according to the 
DIAG36 generator, and avoided possible problems with the soft part of the spectrum 
due to trigger or muon tagging inefficiency. 

• Finally, double-tagged events were rejected by requiring there to be no energy deposit 
exceeding Q.?,xEbeam in the detector arm (defined as ^ = 0°-90° and 6 = 90°-180°) 
opposite that containing the tagged electron. 



Using the high redundancy of the trigger pO|] , the trigger inefficiency was found to be 
negligible for these events. 
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In order to improve the measurements of the tagged electron parameters (energy and 
angles), the following procedures were used. 

1. To avoid edge effects, the tagged electron was required to lie in the polar angle range 
3° <9 < 7.6° (172.4° <9 < 177°) for the SAT, 2.5° < ^ < 9° (171° <9 < 177.5°) 
for the STIC, or 11° < ^ < 35° (145° <0 < 169°) for the FEMC. 

2. To improve the 6 measurements in the SAT, which had a limited granularity, the 
radial position of the cluster was corrected using the function found from the com- 
parison of the experimental radial distribution for Bhabha events with the theoret- 
ical one based on a 1/9^ cross section dependence (Fig. ^). This improved the 
resolution from 6.0% to 2.9%. 

3. To improve the 6 measurements in the SAT and STIC, their alignments were checked 
using Bhabha event samples. The detector on the electron side had a mask in front 
of it to better define the acceptance at low 6. From the number of Bhabha events as a 
function of the electron azimuthal angle 0i, it was possible to find the displacement of 
the mask relative to the beam line. The alignment on the opposite side was checked 
by looking at the difference of the measured polar angles 9tag — Ountag of the scattered 
electron and positron as a function of the positron azimuthal angle 02 (Fig- §)• The 
dependencies observed were used to correct the measured polar angles. The errors 
of the fitted parameters were taken as uncertainties of the procedure, contributing 
0.5% uncertainty on low values of Q^. 

4. A more accurate value of the tagged electron energy Efag was calculated from the 
requirements of energy and longitudinal momentum conservation in the event: 

J-, -P/U/^ cos -|- (2il/beam ^t^l^ ^'^^^untag /-,% 

COS Ountag " COS Otag 

where P^^, E^^ and 0^^ are the momentum, energy and polar angle of the muon 
system, and Ountag is the polar angle of the untagged electron, assumed to be 
or TT. The improvement due to this method can be seen in Fig. |^, obtained from 
simulation, where the difference between the reconstructed and true (generated) tag 
energy Etag — Ef^g is shown as a function of the tag angle Otag using both the direct 
measurement of Etag and this method. 



6 Background 

The following sources of background to the ^~ event samples were considered: 

• e^e~ — s> e^e~r+r^ with a r decay product identified as a muon. The background 
from this process was found to be (1.2±0.2)% for the SAT and STIC tagged samples 
and (5.7±1.1)% for the FEMC, where the errors quoted are statistical. 

• e"'"e'^ r+r^(7) with a hard radiated photon or a r decay product faking a tagged 
electron. This background was found to be negligible for the SAT and STIC samples, 
and (8.9±1.9)% for the FEMC, after taking into account the on-peak versus off-peak 
luminosity distribution of the data. 

• e"'"e~ — >• /x+/x~(7) with the radiated photon faking a tagged electron. This was found 
to be negligible due to the 30 GeV cut on the sum of the muon momenta. 

• e^e~ e+e^TT+TT^ with a pion misidentified as a muon. The ratio of the cross 
sections for pion pair and muon pair production in two-photon interactions falls to 
(1-5)% if the invariant mass of the produced pair is above 2.0 GeV/c^ |26|. With the 
muon identification criteria described above, the probability to misidentify a pion 
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as a muon was below 1.5% (depending on the pion momentum), so this background 
was also negligible for all samples. 

• other e^e^ e+e^ + hadrons processes. These were also found to be negligible for 
all event samples. 

• untagged e^e~ —>■ e^e^fi^fi^ in coincidence with an off-momentum electron faking 
a tagged electron. The off-momentum electrons are beam electrons that have scat- 
tered off residual gas molecules inside the beam pipe. Using a method similar to 
the one described in ^7j, this background was estimated from Z° /i/x events in 
coincidence with a similar off-momentum electron, multiplying by the ratio of the 
dimuon production cross sections from untagged two-photon interactions and from 
Z° decays, and was also found to be negligible. 



7 Results 

The numbers of selected data events after background subtraction are compared with 
the predictions of the signal Monte Carlo simulations in Table 1. The ranges shown are 
calculated given the angular coverage of the detectors and the cut on the tag energy, and 
the average values < > are taken from the data. Figs. |^-^ present the distributions of 
a standard set of observables for events tagged by the SAT, STIC and FEMC respectively. 



Tagging detector 
g2 range (GeVVc^) 
< g2 >*) (GeVVc^) 


SAT 
3.4-36.6 
13.0 


STIC 
2.4-51.2 
12.1 


FEMC 
45.9-752.8 
120.0**) 


data 

BDKRC simulation 
DIAG36 simulation 


1357±37 
1362±14 
1298±25 


2875±54 
2884±22 
2785±55 


239±18 
250±6 
236±13 



*) After requiring Etag > 0.75 x Ebeam (see text) 

**) For events with Otag < 25° {9tag > 155°) (see text) 



Table 1: Numbers of selected events after background subtraction. 

Table 1 and Figs. |§ - 1§ show that the BDKRC and DIAG36 generators produce similar 
kinematical distributions, but DIAG36 gives somewhat lower numbers of selected events. 
In the kinematical region under study, the contribution of the additional diagrams in 
DIAG36 was found to be very small (see section 5). This difference (if real) should 
therefore be attributed to the effect of radiative corrections. The BDKRC generator was 
therefore used for the structure function studies below. 
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7.1 Extraction of K 

To extract F2 , the experimental x distribution was divided by the Monte Carlo dis- 
tribution weighted by the factor a/F2{x,Q'^), where F2{x,Q'^) can be obtained from a 
simulated event sample using either a generator producing events according to a given F2 



or the photon flux approach described, for example, in and briefly outlined below. 
It follows from Eqs. (1-3) that, neglecting the small contribution from terms: 
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where the weight VV(x, Q^) is given by 



2 'untag 



untag 



To calculate the integration limits, the fractional energy of the target photon is extracted 
from the expressions for x and W^^: 

_ 2Q'^/sx - 2Q'^/s + cos 9 + cosg^n^g - 0:^3(1 + cos 9) 

COS 9 + COS Ountag + Xtag{l - COS 9) 

where 

COS 9 = sin 9tag sin 9untag COS(A0) - COS 9tag COS 9untag, (H) 

and A0 is the azimuthal angle between the scattered e"*" and e~. In the single tag 
approximation, 9untag — so that (10) becomes: 

_ 2Q'^/SX - 2QVs - cos + 1 - Xtagjl - cos 9tag) 

- COS 9tag + 1 + Xtag{l + COS 9tag) 

The maximum and minimum Xuntag values correspond to the minimum and maximum 
Xtag values, and these result from the tagging conditions: 

xTr = mm\l-^^A r- , x"!r = max^.^ ^ ' 



^^"^^ I E,^aJ .sin2(^Sr/2)J' I ^ ' ss.n\9Z^/2)\ 

(13) 

where -E^™ is the lower cut on the tag energy and 9^^{9'^g^) is the lower (upper) angular 
acceptance of the tagging device. -E^^" was increased from 0.6 x Eheam to 0.75 x Ebeam 
in order to keep the contribution small. 

Fig. 1^ shows the F2 [x) values obtained by both methods for a simulated event sample 
with STIC tagging conditions, demonstrating that they give similar results. 

A fit to the QED prediction 1T|J2§ 



= I {x^ + (1 - xf) In , ^^J^^ - 1 + 8x(l - x) ^'^^^ 



ml + P^x(l — x) ml + P'^xil — x) 



where terms of order m^/Q^ are neglected gives values of the effective average target pho- 
ton virtuality of 0.022±0.007 and 0.026±0.006 GeV^ for the first and second methods 
respectively, the errors quoted being statistical. For the SAT tagged events the first 
method, which was chosen for the further analysis, gives P^=0.032±0.007 GeV^, demon- 
strating the need to take the target photon virtuality into account in studies of F2 . 

The extracted structure function < F2{x,Q^) >, transformed to F2{x,< >) us- 
ing the ratio F2{x,< Q"^ >)/ < F2{x,Q^) > predicted by QED, is shown in Table 2 
and Fig. |IU|, which present the weighted combination of the SAT and STIC results with 
< >= 12.5 GeVVc^ and the FEMC result with < >= 120 GeVVc^ The FEMC 
sample included only events with 9tag below 25° (above 155°) in order to exclude the re- 
gion with large background contamination (Fig. |^), and the contribution from diagrams 
other than the multiperipheral one predicted by the BDKRC generator was subtracted. 
The structure function values have been corrected to the centres of the x bins by multi- 
plying the measured average values of F2 for each x bin by the ratio of the value of F2 in 
the centre of the bin to the its average value over the bin predicted by QED. Systematic 
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errors due to the resolutions in and x have been evaluated in simulation by varying 
these variables according to their resolutions and checking the effect on F^. The role of 
the observed discrepancy between the data and simulation in some Qtag intervals (Fig. ||) 
was checked by weighting the contributions of events in those intervals according to their 
Qtag values when producing the x distribution. The largest contribution to the systematic 
error comes from the resolution. 

Fits to the QED prediction (14) give = o.025 ± 0.005 and 0.073 ± 0.056 GeV^ for 
the samples with low and high respectively, in good agreement with the Monte Carlo 
prediction. 





<0.1 


0.1-0.2 


0.2-0.3 


0.3-0.4 


0.4-0.5 


0.5-0.6 


0.6-0.7 


0.7-0.8 


>0.8 


stat. error 
syst. error 


0.106 
±0.008 
±0.023 


0.273 
±0.012 
±0.012 


0.426 
±0.017 
±0.012 


0.515 
±0.021 
±0.012 


0.573 
±0.024 
±0.004 


0.645 
±0.029 
±0.003 


0.743 
±0.038 
±0.021 


0.942 
±0.060 
±0.053 


1.152 
±0.112 
±0.094 





<0.2 


0.2-0.4 


0.4-0.6 


0.6-0.8 


>0.8 


stat. error 
syst. error 


0.387 
±0.214 
±0.015 


0.464 
±0.133 
±0.051 


0.673 
±0.138 
±0.049 


0.984 
±0.162 
±0.026 


1.508 
±0.231 
±0.044 



Table 2: The measured structure function for < >= 12.5 (upper table) and 
120 GeVVc^ (lower table). 



7.2 Azimuthal correlations 

In order to increase the observed azimuthal correlations of the final state particles, 
only events with 20° < 9* < 160° have been considered. Taking into account the anti- 
symmetry of Fj in cos 6'*, events with cos^^* <0 and cos^^* >0 have been combined using 
the transformation X — ^ vr — 

The selected samples have been corrected for detector acceptance and efficiency us- 
ing either bin- by-bin corrections over a two-dimensional grid of x ^i-nd 6'*, or a three- 



dimensional unfolding p9[ in the space of the variables X) ^* cind x. The corrected 



distributions (Fig. [11]) were fitted to the expression: 

dN/dx = C (1 ± Pi cos X ± P2 cos 2x) (15) 

where Pi and P2 are closely related to F^jF^ and F'^jFl-, c.f. Eq. (6). The combined 
results were obtained by refitting the weighted sums of corrected distributions for the 



SAT and STIC samples (Fig. |T2D. The parameters determined from the fit are shown in 
Table 3. 

The systematic effects were estimated using simulated events, varying the variables 
Q^, H^77, x^ 6* and x according to their resolution, and adding the resulting variations of 
the fitted parameters in quadrature. This gave errors on the fitted parameters of about 
0.02. The difference between the results obtained with the two different correction meth- 
ods gave an additional systematic error of 0.02—0.06. 

The results obtained were extrapolated to the full 6* and W^^ ranges using the theo- 
retical correction factors Cp^ and Cp^ shown in Table 3, which were obtained as ratios of 



the QED predicted structure functions [0 calculated for event samples generated in the 



8 



range of 2.4-51.2 GeV^ without and with the selection cuts. The results thus obtained 



for F'\jF2 and \F'^/F2 are shown in Table 3 and Fig. 13. They are in agreement with 



the theoretical predictions |]18| and with the results of other LEP experiments |T3| , P 
(note the factor -1/2 difference of with [O] due to its different definition). 



8 Conclusions 

Muon pair production in single-tagged 77 colhsions has been studied at y/s ~91 GeV 
using data collected by the DELPHI detector at LEP during the years 1992-95. Distri- 
butions of different event variables for ranging from ~2.5 to ~750 GeV^/c'^ are well 
reproduced by a Monte Carlo simulation based on QED. 

The leptonic structure function F2 has been measured for two regions of momentum 
transfer with average values of 12.5 and 120 GeV^/c^. 

Azimuthal correlations of final state particles have also been studied, giving infor- 
mation on additional structure functions F]^ and Fg. The measured ratios F\jF2 and 
Pb/P"! significantly different from zero and consistent with QED expectations. 
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x-interval 


Pi 


Cp, 


SAT 


STIC 


Combined 


X < 0.2 

0.2 < a; < 0.4 
0.4 < a; < 0.6 
X > 0.6 
all X 


0.19 ±0.14 ±0.03 
0.22 ±0.09 ±0.03 
0.13 ±0.09 ±0.05 
-0.41 ±0.10 ±0.07 
-0.03 ±0.05 ±0.03 


0.28 ±0.08 ±0.04 
0.20 ±0.06 ±0.02 
0.02 ±0.07 ±0.05 
-0.26 ±0.07 ±0.05 
-0.03 ±0.03 ±0.03 


0.25±0.08 
0.20±0.05 
0.06±0.07 
-0.31±0.07 
-0.03±0.04 


0.541 

0.701 
0.625 
0.849 
0.605 



x-interval 


P2 


CP2 


SAT 


STIC 


Combined 


X < 0.2 
0.2 <x < 0.4 

0.4 < a; < 0.6 
X > 0.6 
all X 


0.06 ±0.12 ±0.03 
0.13 ±0.08 ±0.03 
0.15 ±0.08 ±0.04 
0.20 ±0.09 ±0.06 
0.13 ±0.05 ±0.02 


-0.01 ±0.08 ±0.03 
0.16 ±0.06 ±0.02 
0.19 ±0.06 ±0.03 
0.30 ±0.06 ±0.04 
0.15 ±0.03 ±0.02 


0.01±0.07 
0.15±0.05 
0.17±0.06 
0.27±0.06 
0.15±0.03 


0.391 
0.512 
0.581 
0.673 
0.570 



x-interval 


Fa/P2 




X < 0.2 
0.2 < a; < 0.4 
OA<x < 0.6 

X > 0.6 
all X 


0.135 ±0.043 
0.140 ±0.035 
0.038 ± 0.044 
-0.263 ±0.059 
-0.018 ±0.024 


0.004 ±0.027 
0.077 ±0.026 
0.099 ±0.035 
0.182 ±0.040 
0.086 ±0.017 



Table 3: Parameters Pi and P2 of the fit to the azimuthal angle distributions for the 
SAT-tagged, STIC-tagged, and combined event samples with = 2.4-51.2 GeV^. The 
first error is statistical and the second is systematic. Cp^ and Cp^ are the correction factors 
to extrapolate the parameters to the full 6* range (see text). The values extracted for 
Fa/ F2 and \F^/ F2 are shown with statistical and systematic errors added in quadrature. 
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Figure 1: The dominant multiperiph- 
eral diagram for the reaction e~^e~ — > 
e^e~ii+IJ,~. Etag and 9 tag are the energy 
and scattering angle of the tagged electron 
or positron. 



Figure 2: Definitions of the angles x 
and 0* in the 77* centre-of-mass sys- 
tem. 
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Figure 3: Distributions of the radial positions of Bhabha electrons in the SAT: left - 
before correcting the radial cluster positions, right - after correction. The line shows the 
result of the fit to the theoretical prediction. 
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Figure 4: Left - number of detected Bhabha events as a function of the azimuthal angle 
for the SAT module with the mask. Right - difference of the polar angles of Bhabha 
electrons measured by the SAT modules as a function of the azimuthal angle. The lines 
show the results of the fits. 
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Figure 5: Difference of reconstructed and true (generated) tag energy Etag versus the tag 
angle 6 tag in simulated events: left - measured value of Etag, right - Etag value calculated 
from the kinematics of the event. 
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DELPHI 






Figure 6: Distributions for the SAT single tagged events: a) Etag/Ebeam, b) 9tag {l8Q°—9tag 
for positrons), c) squared momentum transfer Q^, d) invariant mass of muon pair, e) sum 
of the transverse momenta of the muons, f) value of x. The points correspond to the 
background subtracted data, the solid line to the BDKRC simulation, and the dashed 
line to the DIAG36 simulation. 
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Figure 7: The same as Fig. ^ for the STIC single tagged events. 
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Figure 8: The same as Figs. |] and |^ for the FEMC single tagged events except that the 
background, relatively much larger here than in Figs. |^ and |^ and shown here by the 
hatched histograms, has in this case been added to the simulated distributions rather 
than subtracted from the data. 
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Figure 9: -F7(a;) for < Q"^ >= 12.5 GeV'^/c^ extracted from a simulated STIC tagged 
event sample (points) and the fit to the QED expression (hnes): full circles and sohd 
lines are for the extraction using the simulation with known F2, open circles and dashed 
lines for the extraction using the photon flux approach. 
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Figure 10: ^2 (a;) extracted from the combined SAT and STIC data (< >= 12.5 
GeVVc^ full circles), and from the FEMC data (< >= 120 GeVVc^ open circles). 
Statistical and systematic errors are added in quadrature. The sohd and dashed lines 
show the QED predictions with P^=0.026 GeV^/c^ for the low and high Q"^ samples, 
respectively. 
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Figure 11: Azimuthal angle distributions corrected for detector inefficiency: a) x <0.2, 
b) 0.2< X <0.4, c) 0.4< X <0.6, d) x >0.6, e) all x. The lines correspond to the results of 
the fit. Pull circles and solid line are for the SAT single tagged events, open circles and 
dashed hne for the STIC events. 
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Figure 12: The same as Fig. 



11 



for the combined SAT and STIC samples. 



20 



DELPHI 




Figure 13: Ratios of leptonic structure functions F]^/ F2 (left) and \F'^/F2 (right) aver- 



aged in the 
predictions from 
is equal to its mean value over the x bin. 



range from 2.4 to 51.2 GeV^ as functions of x. The lines show the QED 
The points are plotted at the x values where the QED prediction 



